Deltamethrin and transfluthrin select for distinct transcriptomic responses in the malaria vector Anopheles gambiae.

BACKGROUND: The widespread use of pyrethroid insecticides in Africa has led to the development of strong resistance in Anopheles mosquitoes. Introducing new active ingredients can contribute to overcome this phenomenon and ensure the effectiveness of vector control strategies. Transfluthrin is a polyfluorinated pyrethroid whose structural conformation was thought to prevent its metabolism by cytochrome P450 monooxygenases in malaria vectors, thus representing a potential alternative for managing P450-mediated resistance occurring in the field. In this study, a controlled selection was used to compare the dynamics of resistance between transfluthrin and the widely used pyrethroid deltamethrin in the mosquito Anopheles gambiae. Then, the associated molecular mechanisms were investigated using target-site mutation genotyping and RNA-seq. METHODS: A field-derived line of An. gambiae carrying resistance alleles at low frequencies was used as starting material for a controlled selection experiment. Adult females were selected across 33 generations with deltamethrin or transfluthrin, resulting in three distinct lines: the Delta-R line (selected with deltamethrin), the Transflu-R line (selected with transfluthrin) and the Tiassale-S line (maintained without selection). Deltamethrin and transfluthrin resistance levels were monitored in each selected line throughout the selection process, as well as the frequency of the L1014F kdr mutation. At generation 17, cross-resistance to other public health insecticides was investigated and transcriptomes were sequenced to compare gene transcription variations and polymorphisms associated with adaptation to each insecticide. RESULTS: A rapid increase in resistance to deltamethrin and transfluthrin was observed throughout the selection process in each selected line in association with an increased frequency of the L1014F kdr mutation. Transcriptomic data support a broader response to transfluthrin selection as compared to deltamethrin selection. For instance, multiple detoxification enzymes and cuticle proteins were specifically over-transcribed in the Transflu-R line including the known pyrethroid metabolizers CYP6M2, CYP9K1 and CYP6AA1 together with other genes previously associated with resistance in An. gambiae. CONCLUSION: This study confirms that recurrent exposure of adult mosquitoes to pyrethroids in a public health context can rapidly select for various resistance mechanisms. In particular, it indicates that in addition to target site mutations, the polyfluorinated pyrethroid transfluthrin can select for a broad metabolic response, which includes some P450s previously associated to resistance to classical pyrethroids. This unexpected finding highlights the need for an in-depth study on the adaptive response of mosquitoes to newly introduced active ingredients in order to effectively guide and support decision-making programmes in malaria control.

De novo determination of mosquitocidal Cry11Aa and Cry11Ba structures from naturally-occurring nanocrystals.

Cry11Aa and Cry11Ba are the two most potent toxins produced by mosquitocidal Bacillus thuringiensis subsp. israelensis and jegathesan, respectively. The toxins naturally crystallize within the host; however, the crystals are too small for structure determination at synchrotron sources. Therefore, we applied serial femtosecond crystallography at X-ray free electron lasers to in vivo-grown nanocrystals of these toxins. The structure of Cry11Aa was determined de novo using the single-wavelength anomalous dispersion method, which in turn enabled the determination of the Cry11Ba structure by molecular replacement. The two structures reveal a new pattern for in vivo crystallization of Cry toxins, whereby each of their three domains packs with a symmetrically identical domain, and a cleavable crystal packing motif is located within the protoxin rather than at the termini. The diversity of in vivo crystallization patterns suggests explanations for their varied levels of toxicity and rational approaches to improve these toxins for mosquito control.

Genomic Shifts, Phenotypic Clines, and Fitness Costs Associated With Cold Tolerance in the Asian Tiger Mosquito.

Climatic variation is a key driver of genetic differentiation and phenotypic traits evolution, and local adaptation to temperature is expected in widespread species. We investigated phenotypic and genomic changes in the native range of the Asian tiger mosquito, Aedes albopictus. We first refine the phylogeographic structure based on genome-wide regions (1,901 double-digest restriction-site associated DNA single nucleotide polymophisms [ddRAD SNPs]) from 41 populations. We then explore the patterns of cold adaptation using phenotypic traits measured in common garden (wing size and cold tolerance) and genotype-temperature associations at targeted candidate regions (51,706 exon-capture SNPs) from nine populations. We confirm the existence of three evolutionary lineages including clades A (Malaysia, Thailand, Cambodia, and Laos), B (China and Okinawa), and C (South Korea and Japan). We identified temperature-associated differentiation in 15 out of 221 candidate regions but none in ddRAD regions, supporting the role of directional selection in detected genes. These include genes involved in lipid metabolism and a circadian clock gene. Most outlier SNPs are differently fixed between clades A and C, whereas clade B has an intermediate pattern. Females are larger at higher latitude yet produce no more eggs, which might favor the storage of energetic reserves in colder climate. Nondiapausing eggs from temperate populations survive better to cold exposure than those from tropical populations, suggesting they are protected from freezing damages but this cold tolerance has a fitness cost in terms of egg viability. Altogether, our results provide strong evidence for the thermal adaptation of A. albopictus across its wide temperature range.

Experimental evolution supports the potential of neonicotinoid-pyrethroid combination for managing insecticide resistance in malaria vectors.

The introduction of neonicotinoids for managing insecticide resistance in mosquitoes is of high interest as they interact with a biochemical target not previously used in public health. In this concern, Bayer developed a combination of the neonicotinoid clothianidin and the pyrethroid deltamethrin (brand name Fludora Fusion) as a new vector control tool. Although this combination proved to be efficient against pyrethroid-resistant mosquitoes, its ability to prevent the selection of pyrethroid and neonicotinoid resistance alleles was not investigated. In this context, the objective of this work was to study the dynamics and the molecular mechanisms of resistance of An. gambiae to the separated or combined components of this combination. A field-derived An. gambiae line carrying resistance alleles to multiple insecticides at low frequencies was used as a starting for 33 successive generations of controlled selection. Resistance levels to each insecticide and target site mutation frequencies were monitored throughout the selection process. Cross resistance to other public health insecticides were also investigated. RNA-seq was used to compare gene transcription variations and polymorphisms across all lines. This study confirmed the potential of this insecticide combination to impair the selection of resistance as compared to its two separated components. Deltamethrin selection led to the rapid enrichment of the kdr L1014F target-site mutation. Clothianidin selection led to the over-transcription of multiple cytochrome P450s including some showing high homology with those conferring neonicotinoid resistance in other insects. A strong selection signature associated with clothianidin selection was also observed on a P450 gene cluster previously associated with resistance. Within this cluster, the gene CYP6M1 showed the highest selection signature together with a transcription profile supporting a role in clothianidin resistance. Modelling the impact of point mutations selected by clothianidin on CYP6M1 protein structure showed that selection retained a protein variant with a modified active site potentially enhancing clothianidin metabolism. In the context of the recent deployment of neonicotinoids for mosquito control and their frequent usage in agriculture, the present study highlights the benefit of combining them with other insecticides for preventing the selection of resistance and sustaining vector control activities.

A genomic amplification affecting a carboxylesterase gene cluster confers organophosphate resistance in the mosquito Aedes aegypti: From genomic characterization to high-throughput field detection.

By altering gene expression and creating paralogs, genomic amplifications represent a key component of short-term adaptive processes. In insects, the use of insecticides can select gene amplifications causing an increased expression of detoxification enzymes, supporting the usefulness of these DNA markers for monitoring the dynamics of resistance alleles in the field. In this context, the present study aims to characterize a genomic amplification event associated with resistance to organophosphate insecticides in the mosquito Aedes aegypti and to develop a molecular assay to monitor the associated resistance alleles in the field. An experimental evolution experiment using a composite population from Laos supported the association between the over-transcription of multiple contiguous carboxylesterase genes on chromosome 2 and resistance to multiple organophosphate insecticides. Combining whole genome sequencing and qPCR on specific genes confirmed the presence of a ~100-Kb amplification spanning at least five carboxylesterase genes at this locus with the co-existence of multiple structural duplication haplotypes. Field data confirmed their circulation in South-East Asia and revealed high copy number polymorphism among and within populations suggesting a trade-off between this resistance mechanism and associated fitness costs. A dual-color multiplex TaqMan assay allowing the rapid detection and copy number quantification of this amplification event in Ae. aegypti was developed and validated on field populations. The routine use of this novel assay will improve the tracking of resistance alleles in this major arbovirus vector.

Serial femtosecond crystallography on in vivo-grown crystals drives elucidation of mosquitocidal Cyt1Aa bioactivation cascade.

Cyt1Aa is the one of four crystalline protoxins produced by mosquitocidal bacterium Bacillus thuringiensis israelensis (Bti) that has been shown to delay the evolution of insect resistance in the field. Limiting our understanding of Bti efficacy and the path to improved toxicity and spectrum has been ignorance of how Cyt1Aa crystallizes in vivo and of its mechanism of toxicity. Here, we use serial femtosecond crystallography to determine the Cyt1Aa protoxin structure from sub-micron-sized crystals produced in Bti. Structures determined under various pH/redox conditions illuminate the role played by previously uncharacterized disulfide-bridge and domain-swapped interfaces from crystal formation in Bti to dissolution in the larval mosquito midgut. Biochemical, toxicological and biophysical methods enable the deconvolution of key steps in the Cyt1Aa bioactivation cascade. We additionally show that the size, shape, production yield, pH sensitivity and toxicity of Cyt1Aa crystals grown in Bti can be controlled by single atom substitution.

Combining genetic crosses and pool targeted DNA-seq for untangling genomic variations associated with resistance to multiple insecticides in the mosquito Aedes aegypti.

In addition to combating vector-borne diseases, studying the adaptation of mosquitoes to insecticides provides a remarkable example of evolution-in-action driving the selection of complex phenotypes. Actually, most resistant mosquito populations show multi-resistance phenotypes as a consequence of the variety of insecticides employed and of the complexity of selected resistance mechanisms. Such complexity makes the identification of alleles conferring resistance to specific insecticides challenging and prevents the development of molecular assays to track them in the field. Here we showed that combining simple genetic crosses with pool targeted DNA-seq can enhance the specificity of resistance allele's detection while maintaining experimental work and sequencing effort at reasonable levels. A multi-resistant population of the mosquito Aedes aegypti was exposed to three distinct insecticides (deltamethrin, bendiocarb and fenitrothion), and survivors to each insecticide were crossed with a susceptible strain to generate three distinct lines. F2 individuals from each line were then segregated based on their survival to two insecticide doses. Hundreds of genes covering all detoxifying enzymes and insecticide targets together with more than 7,000 intergenic regions equally spread over mosquito genome were sequenced from pools of F0 and F2 individuals unexposed or surviving insecticide. Differential coverage analysis identified 39 detoxification enzymes showing an increased gene copy number in association with resistance. Combining an allele frequency filtering approach with a Bayesian F ST-based genome scan identified multiple genomic regions showing strong selection signatures together with 50 nonsynonymous variations associated with resistance. This study provides a simple and cost-effective approach to improve the specificity of resistance allele's detection in multi-resistant populations while reducing false positives frequently arising when comparing populations showing divergent genetic backgrounds. The identification of novel DNA resistance markers opens new opportunities for improving the tracking of insecticide resistance in the field.

Predicting the success of an invader: Niche shift versus niche conservatism.

Invasive species can encounter environments different from their source populations, which may trigger rapid adaptive changes after introduction (niche shift hypothesis). To test this hypothesis, we investigated whether postintroduction evolution is correlated with contrasting environmental conditions between the European invasive and source ranges in the Asian tiger mosquito Aedes albopictus. The comparison of environmental niches occupied in European and source population ranges revealed more than 96% overlap between invasive and source niches, supporting niche conservatism. However, we found evidence for postintroduction genetic evolution by reanalyzing a published ddRADseq genomic dataset from 90 European invasive populations using genotype-environment association (GEA) methods and generalized dissimilarity modeling (GDM). Three loci, among which a putative heat-shock protein, exhibited significant allelic turnover along the gradient of winter precipitation that could be associated with ongoing range expansion. Wing morphometric traits weakly correlated with environmental gradients within Europe, but wing size differed between invasive and source populations located in different climatic areas. Niche similarities between source and invasive ranges might have facilitated the establishment of populations. Nonetheless, we found evidence for environmental-induced adaptive changes after introduction. The ability to rapidly evolve observed in invasive populations (genetic shift) together with a large proportion of unfilled potential suitable areas (80%) pave the way to further spread of Ae. albopictus in Europe.

Bacterial microbiota of Aedes aegypti mosquito larvae is altered by intoxication with Bacillus thuringiensis israelensis.

BACKGROUND: Insect microbiota is a dynamic microbial community that can actively participate in defense against pathogens. Bacillus thuringiensis (Bt) is a natural entomopathogen widely used as a bioinsecticide for pest control. Although Bt's mode of action has been extensively studied, whether the presence of microbiota is mandatory for Bt to effectively kill the insect is still under debate. An association between a higher tolerance and a modified microbiota was already evidenced but a critical point remained to be solved: is the modified microbiota a cause or a consequence of a higher tolerance to Bt? METHODS: In this study we focused on the mosquito species Aedes aegypti, as no work has been performed on Diptera on this topic to date, and on B. thuringiensis israelensis (Bti), which is used worldwide for mosquito control. To avoid using antibiotics to cure bacterial microbiota, mosquito larvae were exposed to an hourly increasing dose of Bti during 25 hours to separate the most susceptible larvae dying quickly from more tolerant individuals, with longer survival. RESULTS: Denaturing gradient gel electrophoresis (DGGE) fingerprinting revealed that mosquito larval bacterial microbiota was strongly affected by Bti infection after only a few hours of exposure. Bacterial microbiota from the most tolerant larvae showed the lowest diversity but the highest inter-individual differences. The proportion of Bti in the host tissue was reduced in the most tolerant larvae as compared to the most susceptible ones, suggesting an active control of Bti infection by the host. CONCLUSIONS: Here we show that a modified microbiota is associated with a higher tolerance of mosquitoes to Bti, but that it is rather a consequence of Bti infection than the cause of the higher tolerance. This study paves the way to future investigations aiming at unraveling the role of host immunity, inter-species bacterial competition and kinetics of host colonization by Bti that could be at the basis of the phenotype observed in this study.

Receptors are affected by selection with each Bacillus thuringiensis israelensis Cry toxin but not with the full Bti mixture in Aedes aegypti.

Bacillus thuringiensis israelensis (Bti) toxins are increasingly used for mosquito control, but little is known about the precise mode of action of each of these toxins, and how they interact to kill mosquito larvae. By using RNA sequencing, we investigated change in gene transcription level and polymorphism variations associated with resistance to each Bti Cry toxin and to the full Bti toxin mixture in the dengue vector Aedes aegypti. The up-regulation of genes related to chitin metabolism in all selected strain suggests a generalist, non-toxin-specific response to Bti selection in Aedes aegypti. Changes in the transcription level and/or protein sequences of several putative Cry toxin receptors (APNs, ALPs, α-amylases, glucoside hydrolases, ABC transporters) were specific to each Cry toxin. Selective sweeps associated with Cry4Aa resistance were detected in 2 ALP and 1 APN genes. The lack of selection of toxin-specific receptors in the Bti-selected strain supports the hypothesis that Cyt toxin acts as a receptor for Cry toxins in mosquitoes.

Alkaline phosphatases are involved in the response of Aedes aegypti larvae to intoxication with Bacillus thuringiensis subsp. israelensis†Cry toxins.

Bacillus thuringiensis subsp. israelensis (Bti) is a natural pathogen of dipterans widely used as a biological insecticide for mosquito control. To characterize the response of mosquitoes to intoxication with Bti, the transcriptome profile of Bti-exposed susceptible Aedes aegypti larvae was analysed using Illumina RNA-seq. Gene expression of 11 alkaline phosphatases (ALPs) was further investigated by real time quantitative polymerase chain reaction and ALP activity was measured in the susceptible strain and in four strains resistant to a single Bti Cry toxin or to Bti. These strains were unexposed or exposed to their toxin of selection. Although all resistant strains constitutively exhibited a higher level of transcription of ALP genes than the susceptible strain, they showed a lower total ALP activity. The intoxication with different individual Cry toxins triggered a global pattern of ALP gene under-transcription in all the one-toxin-resistant strains but involving different specific sets of ALPs in each resistant phenotype. Most of the ALPs involved are not known Cry-binding proteins. RNA interference experiment demonstrated that reducing ALP expression conferred increased the survival of larvae exposed to Cry4Aa, confirming the involvement of ALP in Cry4Aa toxicity.

Impact of agriculture on the selection of insecticide resistance in the malaria vector Anopheles gambiae: a multigenerational study in controlled conditions.

BACKGROUND: Resistance of mosquitoes to insecticides is mainly attributed to their adaptation to vector control interventions. Although pesticides used in agriculture have been frequently mentioned as an additional force driving the selection of resistance, only a few studies were dedicated to validate this hypothesis and characterise the underlying mechanisms. While insecticide resistance is rising dramatically in Africa, deciphering how agriculture affects resistance is crucial for improving resistance management strategies. In this context, the multigenerational effect of agricultural pollutants on the selection of insecticide resistance was examined in Anopheles gambiae. METHODS: An urban Tanzanian An. gambiae population displaying a low resistance level was used as a parental strain for a selection experiment across 20 generations. At each generation larvae were selected with a mixture containing pesticides and herbicides classically used in agriculture in Africa. The resistance levels of adults to deltamethrin, DDT and bendiocarb were compared between the selected and non-selected strains across the selection process together with the frequency of kdr mutations. A microarray approach was used for pinpointing transcription level variations selected by the agricultural pesticide mixture at the adult stage. RESULTS: A gradual increase of adult resistance to all insecticides was observed across the selection process. The frequency of the L1014S kdr mutation rose from 1.6% to 12.5% after 20 generations of selection. Microarray analysis identified 90 transcripts over-transcribed in the selected strain as compared to the parental and the non-selected strains. Genes encoding cuticle proteins, detoxification enzymes, proteins linked to neurotransmitter activity and transcription regulators were mainly affected. RT-qPCR transcription profiling of candidate genes across multiple generations supported their link with insecticide resistance. CONCLUSIONS: This study confirms the potency of agriculture in selecting for insecticide resistance in malaria vectors. We demonstrated that the recurrent exposure of larvae to agricultural pollutants can select for resistance mechanisms to vector control insecticides at the adult stage. Our data suggest that in addition to selected target-site resistance mutations, agricultural pollutants may also favor cuticle, metabolic and synaptic transmission-based resistance mechanisms. These results emphasize the need for integrated resistance management strategies taking into account agriculture activities.

GeBP/GPL transcription factors regulate a subset of CPR5-dependent processes.

The CONSTITUTIVE EXPRESSOR OF PATHOGENESIS-RELATED GENES5 (CPR5) gene of Arabidopsis (Arabidopsis thaliana) encodes a putative membrane protein of unknown biochemical function and displays highly pleiotropic functions, particularly in pathogen responses, cell proliferation, cell expansion, and cell death. Here, we demonstrate a link between CPR5 and the GLABRA1 ENHANCER BINDING PROTEIN (GeBP) family of transcription factors. We investigated the primary role of the GeBP/GeBP-like (GPL) genes using transcriptomic analysis of the quadruple gebp gpl1,2,3 mutant and one overexpressing line that displays several cpr5-like phenotypes including dwarfism, spontaneous necrotic lesions, and increased pathogen resistance. We found that GeBP/GPLs regulate a set of genes that represents a subset of the CPR5 pathway. This subset includes genes involved in response to stress as well as cell wall metabolism. Analysis of the quintuple gebp gpl1,2,3 cpr5 mutant indicates that GeBP/GPLs are involved in the control of cell expansion in a CPR5-dependent manner but not in the control of cell proliferation. In addition, to our knowledge, we provide the first evidence that the CPR5 protein is localized in the nucleus of plant cells and that a truncated version of the protein with no transmembrane domain can trigger cpr5-like processes when fused to the VP16 constitutive transcriptional activation domain. Our results provide clues on how CPR5 and GeBP/GPLs play opposite roles in the control of cell expansion and suggest that the CPR5 protein is involved in transcription.

Early stages of Golgi vesicle and tubule formation require diacylglycerol.

We have investigated the role for diacylglycerol (DAG) in membrane bud formation in the Golgi apparatus. Addition of propranolol to specifically inhibit phosphatidate phosphohydrolase (PAP), an enzyme responsible for converting phosphatidic acid into DAG, effectively prevents formation of membrane buds. The effect of PAP inhibition on Golgi membranes is rapid and occurs within 3 min. Removal of the PAP inhibitor then results in a rapid burst of buds, vesicles, and tubules that peaks within 2 min. The inability to form buds in the presence of propranolol does not appear to be correlated with a loss of ARFGAP1 from Golgi membranes, as knockdown of ARFGAP1 by RNA interference has little or no effect on actual bud formation. Rather, knockdown of ARFGAP1 results in an increase in membrane buds and a decrease of vesicles and tubules suggesting it functions in the late stages of scission. How DAG promotes bud formation is discussed.

GeBP and GeBP-like proteins are noncanonical leucine-zipper transcription factors that regulate cytokinin response in Arabidopsis.

Understanding the role of transcription factors (TFs) is essential in reconstructing developmental regulatory networks. The plant-specific GeBP TF family of Arabidopsis thaliana (Arabidopsis) comprises 21 members, all of unknown function. A subset of four members, the founding member GeBP and GeBP-like proteins (GPL) 1, 2, and 3, shares a conserved C-terminal domain. Here we report that GeBP/GPL genes represent a newly defined class of leucine-zipper (Leu-zipper) TFs and that they play a redundant role in cytokinin hormone pathway regulation. Specifically, we demonstrate using yeast, in vitro, and split-yellow fluorescent protein in planta assays that GeBP/GPL proteins form homo- and heterodimers through a noncanonical Leu-zipper motif located in the C-terminal domain. A triple loss-of-function mutant of the three most closely related genes gebp gpl1 gpl2 shows a reduced sensitivity to exogenous cytokinins in a subset of cytokinin responses such as senescence and growth, whereas root inhibition is not affected. We find that transcript levels of type-A cytokinin response genes, which are involved in the negative feedback regulation of cytokinin signaling, are higher in the triple mutant. Using a GPL version that acts as a constitutive transcriptional activator, we show that the regulation of Arabidopsis response regulators (ARRs) is mediated by at least one additional, as yet unknown, repressor acting genetically downstream in the GeBP/GPL pathway. Our results indicate that GeBP/GPL genes encode a new class of unconventional Leu-zipper TF proteins and suggest that their role in the cytokinin pathway is to antagonize the negative feedback regulation on ARR genes to trigger the cytokinin response.

Carboxyl-terminal disulfide bond of acid sphingomyelinase is critical for its secretion and enzymatic function.

The human acid sphingomyelinase (ASM, EC 3.1.4.12), a lysosomal and secretory protein coded by the sphingomyelin phosphodiesterase 1 (SMPD-1) gene, catalyzes the degradation of sphingomyelin (SM) to ceramide and phosphorylcholine. We examined the structural-functional properties of its carboxyl-terminus (amino acids 462-629), which harbors approximately 1/3 of all mutations discovered in the SMPD-1 gene. We created four naturally occurring mutants (DeltaR608, R496L, G577A, and Y537H) and five serial carboxyl-terminal deletion mutants (N620, N590, N570, N510, and N490). Transient transfection of the His/V5-tagged wild-type and mutant recombinant ASM in Chinese hamster ovary cells showed that all the mutants were normally expressed. Nonetheless, none of them, except the smallest deletion mutant N620 that preserved all post-translational modifications, were found capable of secretion to the medium. Furthermore, only the N620 conserved functional integrity (100% activity of the wild type); all other mutants completely lost the ability to catalyze SM hydrolysis. Importantly, cell surface biotinylation revealed that mutant DeltaR608 transfected CHO cells and fibroblasts from a compound heterozygous Niemann-Pick disease type B (NPD-B) patient (DeltaR608 and R441X) have defective translocation to the plasma membrane. Furthermore, we demonstrated that the DeltaR608 and N590 were trapped in the endoplasmic reticulum (ER) quality control checkpoint in contrast to the wild-type lysosomal localization. Interestingly, while the steady-state levels of ubiquitination were minimal for the wild-type ASM, a significant amount of Lys63-linked polyubiquitinated DeltaR608 and N590 could be purified by S5a-affinity chromatography, indicating an important misfolding in the carboxyl-terminal mutants. Altogether, we provide evidence that the carboxyl-terminus of the ASM is crucial for its protein structure, which in turns dictates the enzymatic function and secretion.

The expression of cell proliferation-related genes in early developing flowers is affected by a fruit load reduction in tomato plants.

Changes in photoassimilate partitioning between source and sink organs significantly affect fruit development and size. In this study, a comparison was made of tomato plants (Solanum lycopersicum L.) grown under a low fruit load (one fruit per truss, L1 plants) and under a standard fruit load (five fruits per truss, L5 plants), at morphological, biochemical, and molecular levels. Fruit load reduction resulted in increased photoassimilate availability in the plant and in increased growth rates in all plant organs analysed (root, stem, leaf, flower, and fruit). Larger flower and fruit size in L1 plants were correlated with higher cell number in the pre-anthesis ovary. This was probably due to the acceleration of the flower growth rate since other flower developmental parameters (schedule and time-course) remained otherwise unaffected. Using RT-PCR, it was shown that the transcript levels of CYCB2;1 (cyclin) and CDKB2;1 (cyclin-dependent kinase), two mitosis-specific genes, strongly increased early in developing flower buds. Remarkably, the transcript abundance of CYCD3;1, a D-type cyclin potentially involved in cell cycle regulation in response to mitogenic signals, also increased by more than 5-fold at very early stages of L1 flower development. By contrast, transcripts from fw2.2, a putative negative regulator of cell division in tomato fruit, strongly decreased in developing flower bud, as confirmed by in situ hybridization studies. Taken together, these results suggest that changes in carbohydrate partitioning could control fruit size through the regulation of cell proliferation-related genes at very early stages of flower development.